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TECHNICAL NOTE 2511 


CALCULATION OF HIGHER APPROXIMATIONS FOR TWO-DIMENSIONAL 
COMPRESSIBLE FLOW BY A SIMPLIFIED ITERATrON PROCESS 
By W. H. Braun and M. M. Klein 


SUMMARY 

The iteration equations for a simplified solution of the nonlinear 
compressihle-flow eq.uation are developed and applied to two profiles, 
the ellipse and the Kaplan section. Limitations imposed on the solu- 
tions hy simplifications in the differential equation and the boundary 
conditions are discussed. The Velocities near the midchord, the ’criti- 
cal and potential limit free-stream Mach numbers, the extent of 
Isentropic supersonic regions are calculated to foxn* approximations 
for the Kaplan section and to six approximations for the ellipse . The 
development of the Iteration equations and the presentation of the 
results are made in conformity with the Karman transonic similarity 
law and comparisons are made with other solutions. 


INTRODUCTION 

Because of its nonlinear natvtre, the partial differential equa- 
tion that describes the flow of a perfect con 5 )re 8 sible fluid has not 
yielded con 5 >letely to any method of solution so far advanced. In 
fact, in such methods as the Iteration procedures and -the variational 
method (reference l), it is not even attempted to find a closed-form 
solution or the general term of a series expansion; rather, a finite 
(and usually small) number of terms of a series expansion are obtained. 
Unfortunately, in such cases, whereas more and more terms are required 
as compressibility effects become more pronounced, the difficulty of 
securing the higher terms becomes inordinately greater. Althou^ the 
difficulties of nonlinearity may be overcome by a transformation to 
the holograph plane and the choice of a convenient gas law (refer- 
ence 2) , a new problem arises in the solution of any boundary- value 
problem because of the distortion of the boundary. 

The solirtions by Iteration are begun by expanding .the velocity 
potential as a series in terms of some flow parameter; it is then pos- 
sible to determine each term from those which precede it by solving a 
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linear equation. 'Wiienever the flow is over a thin hody at relatively 
hi^ speeds, it is preferable to eacpand the potential in terms of some 
parameter characteristic of the hody, such as thickness ratio, camber 
coefficient, or angle of attack. The first term of the solution is a 
constant representing the undisturbed flow and the second is the well- 
known linearized or small perturbation congpressible flow (Prandtl- 
Glauert solution) . The procedure for obtaining successive terms is 
outlined by Ackeret (reference 3) and extended by Kaplan (reference 4), 
iho applies it to a bun® with no stagnation points and to a circular- 
arc profile. The higher-order terms were shown to be of greater con- 
sequence as the Mach number increased. 

Recently Perl and Klein (reference 5) applied the Prandtl- Ackeret 
procedure, without specifying beforehand the parameter of expansion, 
to flows about thin bodies at transonic speeds. In their analysis, 
the flow equation is transformed to a new set of coordinates in which 
it is effectively expanded in powers of 3^ = 1 - Mq^, where Mq is 
the f ree-stream Mach number . In the transonic range ( "*■ O) , the 

p 

lowest order terms in p are dominant, and the transonic part of the 
flow equation is easily extracted. 

The resulting transonic equation does not of Itself preclude size- 
able perturbations as long as fr is chosen sufficiently close to 
zero; however, the Perl-Klein process uses the customary small pertur- 
bation bo\mdary condition, so that the solution of the complete prob- 
lem - differential equation and boundary conditions - is termed the 
transonic limiting solution, limiting in the sense of Mach number near 
unity and geometric parameter near zero . This limiting nature of the 
solution is further emphasized by the form of the expansion parameter, 
which is found to be simply Karman’s transonic similarity parameter 
(reference 6) , a quantity which combines the three definitive character- 
istics of the flow; the ratio of specific heats (that is, the type 
of gas), the f ree-stream Mach number, and the geometric parameter. 

Although the iteration process greatly facilitates the solution 
of a nonlinear problem by substituting for it a series of linear prob- 
lems, it, in turn, raises the difficult question of convergence of the 
series solution- Tbe present analysis iras made at the EACA Lewis lab- 
oratory to extend the calculations necessary to estimate the convergence 
of the iteration process of reference 6 where it is applied to two 
profiles, the ellipse and the Kaplan section- In partlculeir, it trin 
be fo\md that the successive approximations suggest convergence not 
only of completely subsonic flows but also of some Isentropic mixed 
flows . 

There has been some dispute ccmiceming the stability of the second 
type flow (see the resume of these discussions by Sears, reference l); 
one of the latest contributions is a stability investigation by Kuo 
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(reference 8), which indicates that decelerating transonic flows are 
unstable for a certain type disturbance. The result of the instability 
is the formation of a terminal shock in the. supersonic region^ with a 
consequent asymmetiy. Some British measurements (references 9 and lO), 
however, indicate that small supersonic • regions not terminated by shock 
■vraves may exist, or if shocks do occur, it is possible that they may 
'not change- the flow greatly from the isentropic pattern, provided that 
the free-stream Mach number is only sll^tly greater than critical. 

The isentropic mixed flows are therefore dealt with herein under the 
assumption that they are indicative of the natiare of some real flows. 

Another unpleasant feature of iteration methods, in addition to 
the question of convergence, is the large amount of labor involved. In 
obtaining the higher approximations by the present method, it is con- 
venient to make certain simplifications in the differential equation 
so that the labor does not become too great. The justification for 
such changes and the limitations they impose on the solution will be 
discussed as they arise. 


SMALL-PERTUEBATIOW TRMSOMIC FLOW 

In reference 5 it is shorm that the partial differential equation 
governing two-dimensional, irrotational, compressible flow may be writ- 
ten 

(P ‘Pxx 'Pyy^ [3 " ^x '*’y )J 

= [% (2«Px + "Px^) + ^ Mo\^«Pxx + + 2Mo^ (1 + Vxy 

( 1 ) 

•where 
= 1 - 

<p pertiorbatlon potential 
y ratio of specific heats of the gas 
Mn free-stream Mach niunber 


•7 


rM = V (1 +I^Mo2) 
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ani^ the free-stream velocity 1i8ls magnitude unity and the direction of 
positive X. (A complete list of synibols is given in the appendix.) 
If the flow is over an isolated hody \7hoBe surface is defined hy 

y = T g(x) 

where t is the thickness ratio, the boundary conditions on the per- 
turbation potential are 

(py = t( 1 + (pjj) gx(3c) at y = T g(x) 

= <py = 0 at x^ + y^ = » 



For purposes of solution, it is desirable to ^nrite the second-order 
terms of equation (l) in normal form, as may be accomp3j.shed by the 
Prandtl-Glauert transformation. If the transformation is modified 
subtly, however, the coefficient of a nonlinear term of the same 
order in also becomes independent of flow parameters. The modi- 

fied transformation is 

X = X 

n = Py 

F(x,ti) = q>(x,y) 



and it transforms the problem consisting of the differential equation (l) 
and the boundary conditions (2) into 


(l - 2 Fj{.)F2qj+ F^^ 


2Mq" \ 
+ „ FJ^/4- 


20C • ^TfxT] 


2 / 12., 

\ M 

\ M / 




2 p 2 11 

M 


p 2 T] TlTlJ 
M 


( 4 ) 
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= (1 + <Px) Sx at T] = pTg(x) (5) 

Fjj. = F^ = 0 at + Ti2/p2 ^ „ (g) 

2 

Figure 1 is ttie reciprocal of tlie transformation coefficient of 

equation (3) . 


A transition is now made to th,e small -perturbation transonic prob- 
lem. For transonic speeds, the parameter approaches zero so that 

the rigjht side of equation (4) may be neglected. Because of the non- 
linearity of the right side, omission of these terms is especially 
Justified if the perturbations are small. If the body considered is 
thin (l »Tr) and there are no stagnation points, q> 2 j.< <1 and ’the 
velocity <pjj. may be neglected in boundary condition (5) . For flows in 
which stagnation points occvlt, the effect of these stagnation points 
is assumed negligible in the region of interest in the flow field. 
Furthermore, in botmdary condition (5), because the product ^ is a 
very small number and because (as will be shown) the solution for 
equation (v), "vdiich follows, has no singularities on the x-axis, it 
is assmed that the condition (that the velocity normal to the body is 
2 !ero) may be made on the x-axis rather than on the body itself. This 
assvimption is Justified in reference 5. Consequently, the small- 
pertixrbation transonic flow problem is defined by 


(l - 2F^) Fyy + F^^ — 0 


( 7 ) 


■n 3 

TJ = 0 

(8) 

X 

II 

II 

o 

x^ + T)2/p2 = « 

( 9 ) 


For values of F^^ < 
Fjj, > ^ it is hyperbolic, 
turbation velocity, <P^ = 


i, the flow equation (?) is elliptic and for 
2 

In 'figure 2 the parabolic value of the per- 
1 

~ — (or F-y- = — ), is compared -with the excess * 

2Tm ^ 2' 


of the pritlcal velocity over free-stream velocity, that is, the pertur- 
bation velocity at which a flow becomes locally sonic, as calciilated 
from the Bernoulli equation. Fpr such values of Mq for which the two 
curves substantially eigree, the differential equation (?) may be expected 
to describe mixed flows properly. The agreement is acceptable in the 
range 0.9 < Mq < 1.1, Tdiich is approximately the range of accuracy which 
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woiild have heeii. suggested "by equation (4) , It is also ohserved at this 
point that if the error in the houndaiy condition (8) is to he limited 
to 10 percent as con 5 >ared with the exact boundary condition (5), it is 

necessary that <p„ < Thus an examination of both the differential 

-*■ 10 

equation and the boundary condition indicates that flows in the transonic 
domain idLll be described sufficiently accurately only if they lie in the 

interval 0< j 

The floi-rs to be solved herein will be limited to those about bodies 
■I'diich are symmetric about the x- and y-axes. A limitation to subsonic 
free- stream Mach numbers is also necessary in order that the trans- 
formation (8) does not- become imaginary. The results of the present 
method for hi^ subsonic speeds can be applied^ however, to slightly 
supersonic free-stream flows for idiich the detached shock wave is far 
ahead of the body and the Mach number do'VTnstream of the shock is as 
much below unity as the Mach number upstream is above unity. In this 
sense, the analysis is applicable to the entire range of transonic 
velocities . 


ITERATIOir PROCESS 

For convenience in describing the bo^, a transformation is made 
to elliptic coordinates, one of which is cyclic and can indicate posi- 
tion arotind the circumference of the body. These coordinates are com- 
pared ^-Tith cartesian coordinates in figure 3. The equations defining 
the transformation are 




X = cosh s cos t 

( 10 ) 

• 


T) = sinh s sin t 

(3J-) 



p P *? *? 

cosh, s - cos t = sinh s + siji "t 

( 12 ) 


s,ty 


where is the Jacobian or functional determinant of the 

\s,t; 

x,y-coordinates with respect to s and t. The reciprocal of the 
Jacobian of the transformation is plotted in figure 4 as a frmction 
of t for several values of s. 


The first derivatives of the potential in the new system are 
F^j. = (sinh s cos t Fg - cosh s sin t F^)j'"^ 

F^ = (cosh s sin t Fg + sinh s cos t F^)j“^ 


(13) 

(14) 
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whlcli Isecome at the axis of the hody 


F = = — Fa. 

2 C sin t t 


^11 sin t 


(15) 

(16) 


Substitution of these into the hoxmdary condition (8) yields 




F = 
s 




3 


s = 0 


(17) 


Boundaary condition (9) becomes 


^8 = ^t = 0 8 = » (18) 

An exact transformation of the flow equation (7) would result in 
the Inverse Janobian appearing in the nonlinear term, and it would be 
necessary to expand it in a Fourier series before en 5 >loying a method 
of solution involving the separation of variables. The process used 
herein would then req.uire the solution of an infinite number of ordinary 
differential equations. In order to simplify the analysis and to reduce 
•the coii 5 )utational work that would ordinarily accompany such an expan- 
sion, only its first or constant term is retained. That is, in the 
transformation of equation (7) to elliptic coordinates, the Jacobian 
is set equal to 1, which is the first term' in its expansion about the 
point s=0, t=jt/2. Of course, this simplification will restrict 
the region in which the equation is valid. A further restriction on 
the size of the region of validity is made by restricting the calcula- 
tions to velocities on or near the surface of the body; thus s *» 0 in 
the differential equation. 

If the relation 

^xx ^ (^Bs ^tt) 

and equations (13) and (li) for the derivatives of F^- and F.« are 
used, equation (7) with the approximations 

1 

sinh B * 0 


and 


cosh 6 « 1 
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may be in-itten 

^ss + ^tt = - 2 sin^t ^ (sin t F^) 

= - 2 ain^ ^t^tt ~ ^ sin^t cos t F^F^ (19) 

Tliis equation is proposed for use as a description of the flow in a 
region about the origin whose boundaary is set by the criteria 

j-1 - ij ^ 0-10 and s <-0.1. Such a region in the xT)-plane is shoim 
in figure 5j where for comparison a 10-percent -thick (in the xy-plane) 
Kaplan section is included as a typical profile. The length of the 
rectangle is seen to be about 0.3 chord. 

With the change in the region of validity of the dlfferentieil 
equation, there must be a change in the boundary conditions to corres- 
pond to this region. The surface boimdary condition (7) may be 
retained as isj therefore only the outer boundaries (dashed rectangle 
in fig. 5) of the region lack boundary conditions. As a substitute 
for conditions on these boundaries , the boundary condition (18) at 
infinity is taken eis a guide to postulate that, in the restricted 
region, both Fg and F^ must decrease with increasing s . The fav- 
orable comparison with more conrplete solutions shows that this con- 
dition prescribing monotonic character of the potential derivatives is 
sufficient. These comparisons are made in subsequent sections. 

The method of solution of the problem, as formulated by the dif- 
ferential equation (19)., the boundary condition (17), and the condition 
of monotonelty in s, is to write the transformed potential F as an 
infinite series 

12 3 

FHF+F+F + ... ( 20)' 

where it is assumed that the succeeding terms become progressively 
smaller at such a rate as to insure convergence of F and its deriva- 
tives. Definition (20) is placed in the flow equation (19), which is 
then arranged as follows: 
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+ F. 


tt 


fo 


2 2 

+^SS 

3 3 

+^ss + ^tt )= { 


n n 
+^ss ■*■ ^tt 


1 1 


1 1 




- 2 sin^t F^F^^ - 2 sin^t cos t F^F^ 

2 112 12 

- 2 sin^t(F^.F^^ + F^F^^) - 2 sin^t cos t 2 F^^ 


n-1 V n-v 
2 sin\ 2 F^F^^ 
v=l . 


n-1 V n-v 

2 sin^t cos t- 2 F^F^ 

. v=i 


( 21 ) 


This arrangement Indicates that the equation should he divided for pur- 
poses of iteration into a sequence of equations each corresponding to 
one line of equation ( 2l) . In this way all terms idilch have the same 
total superscript value are grouped together; this procedure will later 
be seen equivalent to separating out terms of the same order in the 
Karman similarity parameter. 

The boimdary condition on the body is also divided in a particular 

1 

way among the terms in the series solution. The first term F mates 
the total contribution; 


s " p3 ^ 


s =,0 


( 22 ) 


VThereas the boundary conditions for all other terms are homogeneous c 

n 

Fg = 0 s = 0_ n>2 (23) 

n n 

The condition of monotoneity requires that all F^ and F^ are monotone 
decreasing in s. 


1 

FIRST APPROXIMATION, F 

The first appro xima tion differs from subsequent approximations in 
that it requires the solution of the homogeneous Laplace equation (first 
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line of equation (2l)) "with a nonhomogeneous boundary condition (22). 
Hl^er approximations inTolve the solution of nonbomogeneous , or Poisson^ 
eq[uations ■witb bomogeneous boundary conditions . For this reason and 
to obtain the necessary start in the iteration procedure^ the first 
approximation is treated separately. 

The solution of this first iteration, -which takes account of the 
monotonelty of the derivatives and of t as the cyclic' coordinate, is 

1 “ _ts 

F = Z e (ajj. cos kt + cdjj- sin kt) (24) 

k=l 


Before this form is introduced into the boundary condition (22), the 
function g (0,t) is expanded into a Fourier series and written 

g (0,t) = Z 5^ sin jt (25) 


In conformance -with the restriction to bodies which are symmetric about 
the y-axis £ind about -which there is no circulation, no cosine terms 
appear in the expansion. The index j is limited to odd values 

because of the symmetry of g about t = iJ. The introduction of expres- 

2 

sions (24) and (25) into the boundary condition on the body (22) leads 
to 


- Z k(cjj. cos kt + sin kt) ^ Z cos jt 

k=l J=1 


(26) 


It follows that 


J odd only 




Z 

k=l 


°k ® 


COB kt 


(27) 


k odd only 


1 

The nature of the dependence of F on 



a quan t ity combining 


the entire dependence of the first approximation on the parameters of 
the flow (reference 6), has the consequences immediately observed from 
equation (2l) that 


/» 


2275 



2275 


HA.CA M 2511 


11 



The quantity 





(28) 


J 


is recognized an the similarity parameter of reference 5, and reduces 
to Karman's similarity parameter (reference 6) as 1, Hence, as 

previously claimed, equation (2l) is the ordering in powers of K with 
respect to this parameter. The similarity parameter is plotted as a 
function of Mq for various x in figure 6, -which is reproduced from 
reference 5. 


n n 


Writing F 3 fK?^ gives the perturhation potential, hy the trans- 
formation (3), as 




12 3 

f fK -f- fK^ + : 


■■■) 


- T 

P 


“* -ks 

2 6^6 cos kt 

L\k=l 


^ + f^K^ -1- . . .J 


(29) 


The first term inside the parenthesis represents an incompressible flo-w. 
Hence, for very small values of x, for -which K is also small, the 
transonic potential (29) reduces to the Prandtl-Glauert rule. 


EEGHER APPROXIMATIONS 
n 

Inasmuch as all the F heyond the first satisfy the same general 
form of differential equation and have the same boundary condition, it 
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n 

is possible "to •treat "tlie general form of. F rather than each approxi- 
mation separately. In the follo'wing discussion^ an induction argument 
•win he used -to sho-w that e-very term in the series for the transformed 
potential has the form 


X(p+3) \(p.+3) 

X ' . -3 -3 

F = Z Z 

r=l p=l 


+ ^s + Nj^s^ + . . . )e"^® cos rt 

(30) 

r odd only 


XX X 

•where 1^, and Hip, are constants a nd jj. is -the number of terms 

in the expansion (equation (25)) of q.(t), the function characterizing 
the shape of -the body. The malTi part of -fche argument is to show that 
if equation (30) is true for an X < n, then it is also true for 
X = n. 


From equation (2l), the differential eq.uation satisfied by any 

n 

term F (n > l) is 

n n-1 V n-v n-1 v n-v 

F + F^. = - 2 sin^t 2 F+F^^ - 2 sln2t cos t Z F^-F^ (3l) 

S S 17C ^ y — 


V n-v 

Use of "bhe postulated form (30) for F and F 
Sion for (3l) , 


leads to the expan- 
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''C^+S) (n-v)(^+3) v(ii+3) (n-v)(n+3) 


SB ■*' ^tt 


^ n-1 -3 -3 

- ^ Z Z 2 

V =1 m=l P=1 


-3 

Z 

k=l 


-3 

Z 

r=l 


®^kni "*■ * • * ) "*■ ®^^rp ®^rp "^ • • • ) • 

-(p-Ha}s [T 2 A 2 

e |(-3kr -kr) cos (k+r+l) t + (-3kr +kr) cos (k-r+l)t + 

(3kr^+kr) cos (k-r-l)t + (3kr^-kr) cos (k+r-l)t + 

(kr+kr ) cos (k+r+3)t + (kr-kr^) cos (k+r-3)t + 

(-kr+kr^) cos (k-r+3)t + (-kr-kr^) cos (k-r-3)-^ (32) 


r,n-v 


k, r odd only 


The right side of the foregoing equation is the sum of products involv- 

V V V 

ing the constants %m-> ^km-> * ' *■» Povrers of s, exponentials in 

B, and cosines of angles of multiple t. Accordingly, after all five 
summations have heen accomplished and regrovtpings made according to 
powers of s, exponentials in s, and multiples of t, the equation 
will- he 


n 

F. 


ss 


n 
+ F 


tt 


n(p+3) n(p+3') 
-3- -6 

= 2 2 
h=l j=2 




If a solution 


h odd only 


. ) e COB ht 
(33) 


n(p-+5) 

n -3 n 

F = K^ Z Qji (s) cos ht 

h 


( 34 ) 


h odd only 
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is eiBsumed, tlie coefficient of each, harmonic term in (33) must vanish. 
Thus 


n(n+3) 

2 n gn -6 n n n 

— Oil (s) - h (s) = 2 + sBj^^ + s + 

ds J=1 


. .) e-Js 


(35) 


K> 

K) 




h odd only 

The particular integral of this equation is 


Pt. ~ 2 ^ ^ ~ 




2h 


n 




Be-J" ^ 


,2 ,2 f .2 , 2 a 2 

j -h (o -h ) _ 


n 

“ ®hh® 


-hs/sf 




*) 


n _is 

o/h 


s2 + '43s ^ 63^42h^ 






^ s _s 

^6h ai.2 at-3 


+ -^1 + . . 
4h“ 4h'' 


) 


(36) 


Here the notation 2 refers to a summation over all values of j 


3/11 


-hs 


except 3 = li" complementary function is Gj^^e , , "where Gjj^ may 

he found from the houndary condition on the body ( 23 ) ; 


n 


- h Gjj + 


(— ) " 

\ ds ! a — 0 


n > 2 


(37) 


Then 
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n 



a(^+3) 

-6 n n ^ , _ia 

= Z + sl^j + + 8 + . . .) e 


(38) 


vliere 


= 




n n 


2 j . ^ 6 ^+ 2 h^ 

j 2 .h 2 >*3 (j 2 .i> 2)2 




n 

iiii= 2 




f-J^ 


j 6,1^+211^ , 4.1 




( la n n \ 

hix + %h ”3 + -\ + • • •] 

V 2 h 4 h 4 h / 


n 

1 


II 

•<-:> 

j|n 

■ 

n 

■ 2E 

n 

%i = 

+ BtLh 

n 

1 

n 

II 


■*■ 

11 

/n 1 

n 

Nht = 

" A ®hli 4H 

+ 


4 J 




Cj^j + . . .^ 


4^2 * ^ 4 h 2 


■) 


+ • • • 


4h 


■■) 


n 

°hh = 


■ °lili 61i • • ' 


(39) 



16 


MCA Iff 2511 


(The maxlitiUTn value of h for any u in equations (38) and (39) is 
n(|i+3)-3.) Finally, if equation (38) is placed in the assumed form (34), 
the n^^ approximation becomes 


n(ii+3) 

n - -3 

F = Z 


n(p43) 

-3 

Z 


h=l j=l 




cos ht 


(40) 


h odd only 


■which is exactly the form postulated' in equation (30) . 

V 

It has been sho'wn that if F has the form (30) for v < n, it 
■vrill aiso have "that form for v = n; there remains only to test this 
form on a particular approximation in order to complete the induction 
proof. If -fche second approximation is taken to correspond to the 
n"*"^ approximation of the proposition, a comparison of all pre'vioua 
appr oximat ions (equation (27), in this case) indicates that the n’^^ 
(second) and, hence, all approxlmaticns .subsequent to it must obey 
•bhe f oim (30) . The fact that the method of solution for the first 
appro xima tions differs’ from that of subsequent approximations does 
not Interfere -with the argument. 


The antual mechanj.es of the iteration process have been derived in 
the foregoing development and are outlined by equations (32), (33), 
(39), and (40). For example, after the first approximation has been 
obtained, the dlfferentieil equation for the second may be found by 
setting n = 2 in equation (32) and letting k, r, m, and p assume 

1 

■values in accordance ■wi'fch the form of F. After all the coefficients 
of each harmonic on the right side have been gathered, the differential 
equation -will have the form (33) . The application of the formula (39) 

2 -2 

yields the constants of the solution (40) in terms of the ®hl-> 

11 ' 2 ^ 2 ^ 

. . . of the differential equation. The ^ 

3 *^3 

. . . are then kno-wn and are sufficient to determine the * • 

if the foregoing process is repeated. 


CO 

oi 


MAmroM VELOCin ihcremeht and convergence 


For flo■^^■ over a profile symmetric about the y-sixls, the 
velocity increment occurs at the midchord. At the critical free- stream 
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Mach, number, the flow -will he sonic at the ml dchord, and for hl^er 
free-stream Mach numbers, a local supersonic region •will develop. 

Connected Td-th "the formation ' of this region is> the question of 
convergence of the series obtained from the iteration process to repre- 
sent the maximum velocity increment. This question of convergence is 
especially important because the ma xi mum velocity increment is expected 
to be the first velocity component at any point in the flo'w -to diverge, 
and this divergence signals the falltire of the iteration process to 
serve as a method of solution. If several terms of a series are kno'wn 
and subsequent terms are ansumed to decrease in the same manner, an 
estimate of the domain of convergence may be made as foUo'ws: 

The limiting (t 0, x o) perturbation velocity at the midchord 
is a sum of terms (equations (3), (15), and (30)) 


or 


n 

lim 
xH-o 
y 0 



= 5k^- 


n 


2 

h 


2 h(-l) 


h-1 

2 



lim <Px = 5 2 p 
X -»■ 0 
y -*• 0 


(41) 


(42) 


n • ' n 

The p are seen by definition to depend only on the Lj^j -which intro- 
duce the effect of the shape of the body. The series (42) -will converge 
if the limit of the ratio of successive terms is less than 1 for aH n; 


■ n+1 

- lijn . P < 1 

n ® g^-1 ^ 
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or 

\ 

n-1 

K, = g.l.lD. (43) 

L ° n 

P 

■Hiiere is the ma-xiTTmTii value of the similarity parameter for -which 

convergence of the iteration process -will occxir. 

'TViP so-called potential limit Mach number, or free-stream Mach 
number for -which the solution diverges for any T , is defined hy 





or 




(44) 


SDRFACE 7EL0CITIES 

For small disturbances such that the squares of the perturbation 
velocities may be neglected conpared -vrith unity, the transformed limit- 
ing perturbation velocity on the eUipse is 




n 


K”S Z Lj,j h sin it 
n h j 


or, in -the x,y-plane. 


= ain ' t I 2 2 2 Ljij h sin ht (45) 

n h 3 

Equation (45) -will hold only on that portion of the surface along -which 
J 4 1 (figs . 5 and 6) . 
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The local Mach nimiber at the midchord is given for the small per- 
turhation case hy equation (27) of reference 11: 

m2 - 1 = - + 2Tm q>x (0,0) (46) 


In a previous discussion concerning figure 2 it was seen that the 
hyperhollc region of the limiting equation (7) represents very nearly 
the region of supersonic flows for a perfect fluid. Therefore the 
parahollc curve 



or (p^ 



(47) 


can he used to outline the supersonic region or to find significant 
points on it in the physical plane. 'Whether or not an Itfaration solu- 
tion can correctly describe a supersonic flow may he questioned, for 
as performed here the Iteration procedure is essentially the solution 
of a set of elliptic equations . The applications which follow, althou^ 
not anstrering this question directly, do suggest convergence of the 
solution in the hyperbolic region. 

For any thickness ratio t the critical Mach number, or lowest 
free-stream Mach number at -v-jhich some point in the flow field becomes 
sonic, is given by 






or 


M. 


^ ¥ m/) 


(1 - ^,) 


K 

c 

T 


(48) 


AEPLICATIOW TO EEUPTIC CYLINDER 

The expression that describes the upper sturface of an ellipse 
■whose major axis coincides -^d-th the direction of flow is 

y = T g(x) = T/\/i - x2 

(s ->-0), 


= T sin t 


( 49 ) 
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It follows from equation (27) that 

1 

F = Ke ® cos t (50) 

Subsequent approximations throu^ the sixth have been calculated ^ 

by means of the iteration process previously described. The coefficients n 

n n n n n n n 

means of equations (33) and (39) are presented in table I in matrix 
arrangement. From equation (41), 


"h.1 


and the 


*^hj calculated by 


tp^ (0,0) = I (1 + 0.35237 K + 0.36986 + 0.51836 + 

P 

0.83711 + 1.42014 + . . . ) (5l) 

For comparison, the mR-yiTm-rm disturbance velocity as determined by 
the first three approximations of the iteration solution herein des- 
cribed is plotted (fig. 7) id-th the exact solution of Hantzsche (ref- 
erence 12). Hantzsche *s iteration process is based on an expansion of 
the stream function in powers of the thickness ratio. The agreement 
between the tiro methods appears close enough to justify the calcula- 
tion of hi^er approximations for obtaining a better estimate of the 
potential limit Mach number. 

An estimate is made in figure 8 of the largest value of the 

similarity parameter for which the series for the mR-x-lTTnTm velocity 
Increment ivLll converge. The broken line joins successive estimates 
of obtained by conparing terms by means of a ratio test described 

in the section "MAXIMUM YELOCin IHCEEMEHT AND COMERGENCE." An 
examination of the slopes of the segments of the line indicates that 
the l imi ting ratio of these slopes, for large n, is about l/2. If 
it is assumed that these slopes continue to decrease in a ratio l/2, 
the greatest lower bound of the broken line curve can be calculated 
by a geometric progression. The estimated greatest lower bound of the 
broken line is then about = 0.56. 

That the solution converges even when some of the flow is s^Iper- 
sonlc (hyperbolic) is indicated in figure 9, which shows the contribu- 
tion of each approximation to the mp-x-l rnirm velocity increment for vary- 
ing K. Supersonic flows lie to the right of the parabolic curve and 
the intersection of this curve and the sixth approximation marks the 
value of the parameter K at which the transition from totally sub- 
sonic flows to partly supersonic flows occurs. This transition 



2275 


mCA OM 2511 


21 


point, Kc == 0.395, is also very evident in figure 10 as the intersection 
of the curves, for several values of Mq, of the local Mach number at 
midchord (equations (42) and (46)). 

The surface velocity on the ellipse for three values of K is 
presented in figure 11. The parabolic line, equation (47), is Included 
in each case for comparison. Thus in figure U(a) al l velocities are 
subsonic, -whereas figure H(b) (K = 0.4) sho-ws a partly superspnic 
flow. This is in agreement -tTith the pre-vdous calculation showing 
transition from subsonic flow to supersonic flow at = 0.395. The 
longi-tudinal extent of the symmetric supersonic region is seen to be 
about l/lO chord, based on six approximations. .The 3^t of this series 
of figures (fig. 11(c)) is for a value of K near K^. Although the 
series for the velocity increment may be converging, it is e-vident 
that near the midchord the uncalculated approximations are not neglig- 
ible for this value of K. The -^d-dth of the supersonic region is 
greater than the extent of the region of accuracy of the calculations. 

Curves of critical Mach number and potential limit Mach number 
against thickness ratio have been plotted in figure 12 using eq.ua- 
tlons (48) and (44), respectively. 


APPLICATION TO KAPLAN BIMP 


As a second example, the iteration process is applied to the cal- 
culation of noncirculatory flow past a Kaplan section; this problem 
has been treated in an exact manner to three approximations in refer- ' 
ence 4. The Kaplan section has no stagnation points (fig. 5) so that 
no-^diere is the boundary condition singular. Its parametric repre- 
sentation is (reference 4), 


X = cos t 

y — j (3 sin t - sin 3t) 


’ 


(52) 


Application of equation (27) yields 


1 

P = 



cos t 


-3s 

e 


cos 3t) 


(53) 


The second, third, and fourth approximations have been calculated 
according tb the procedure given in the section "HIGHER APPROXIMATIONS" 

n n 

^d are recorded in table. H as arrays of the constants Lj^j^ nnd 

Njj j . The constants ^h^ represent the nonhomogeneous 

terms of the differential equations (33). 
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The chordvd.se perturhatlon-velocity distrihution at the surface Is 
plotted In figure 13 for comparison -vri-th the result of Kaplan (refer- 
ence 4) for t: = 0.10 and Mq = 0.75. The agreement is prohahly suf- 
ficient from a practical standpoint to Justify the labor-saving limit- 
ing processes that have been made in the iteration process. 

An estimate of the limiting value of the parameter K for -vdiich 
the series for the Tna-idTmTm velocity incr&aent -will converge is made 
in figure 14, "where the ratios of successive coefficients in the series 
for the Tnn-idnnim Velocity increment are plotted for the first four terms . 
The series for (‘Px)max four approximations is calculated from 
equation (4l) to be 

(tp ) = 1 ( J + 0.81762 K + 1.15715 K^ + 2.11564 K^ + . . .) (54) 

' x'max p ^ / 

An analysis similar to that employed on the correspopdlng series for 
the elliptic cylinder indicates that the series (56) -will converge for 
values of K < 0.387. 

The value Kj = 0.387 is considerably above the value of K for 
■vjhich the flow becomes sonic at the midchord, as given in figure 15 by 
the intersection of the fourth approximation curve and the parabolic 
curve. The converging nature of the successive approximations indi- 
cates that this value of = 0.270 "vd-ll be lowered only sli^tly 
by increased accuracy from higher-order terms . 

The variation of the local Mach number corresponding to the maxi- 
mum velocity increment has been calculated by means of equation (46) 
and is plotted against the similarity parameter for several values of 
free-stream Mach number in figure 16. The calculations are based on 
four approximations . The velocity increment on the surface of the bump 
is plotted as a function of distance along the chord in figure 17 . Of 
these, figure 17(a) sho'vrs the only completely subsonic flow. In flg- 
vire 17(b), <K = 0.3 <K^. Hence the flow is everyvdiere convergent, 
but near the midchord there is a small supersonic region represented by 
those values of (Fx)max parabolic curve. For figure 17(c), 

the similarity parameter is chosen slightly larger than the estimated 
ma-ximum for convergence, Kj, as obtained from equation (43) and fig- 
vire 14; and, althou^ successive approximations still decrease monoton- 
Ically, the curves of figure 17(c) should not be expected ultimately 
to converge. 

Figure 18 shows the critical Mach number Mq and the potential 
limit Mach number Mj as functions of thickness ratio. For compari- 
son, and Mj as computed by Kaplan (reference 4) are included. 

The present calculations indicate that isentropic mixed flow is 
restricted to a considerably smaller range of Mach number . 
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COirCLDDIHG HEMAEKS 

An. iteration procedure similar to that en^Jloyed in the Ackeret- 
Prandtl type prqcedure has heen utilized in the calculation of the 
symmetric-type transonic potential flow for an elliptic cylinder and 
a Kaplan section. In order to obtain an estimate of the region of con- 
vergence of this type of- solution, the -\rork has heen carried through 
six approximations for the elliptic cylinder and throu^ four approxi- 
mations for the Kaplan section. 

The resiilts indicate that the iteration solution converges for a 
range of the similarity parameter for -which supersonic regions exist 
in the flo-w field. For these hyperbolic regions, the nonlinear term 
of the potential equation is not, strictly speaking, small compared 
•v^ith either tem of the Laplacian. The iteration procedure presented 
herein may therefore he appll'cahle to a greater range of dis-turhance 
velocities than -would he indicated hy tts me-bhod of formulation, in 
-which the nonlinear term is considered a perturbation on the Laplacian. 

The results have heen presented in conformity "with the transonic 
similarity law. The essential characteristics of a flow may then he 
presented in a simple and direct manner hy classifying the flow -with 
respect -fco significant values of the similarity parameter, for example, 
the critical and potential limit values. 


Lewis Fll^t Propulsion Laboratory 

National Ad-vLsory Committee for Aeronautics 
Cleveland, Ohio, June 26, 1951 
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APPENDIX - SYMBOLS 


The following symbols are used In "this report: 


A,B,C 

c* 

F 

f 

g 

g.l.h. 


constants 
critical velocity 

perturhation potential in xrj-plane 
function of x,ti 

thickness distribution along chord 
greatest lower hound 



Kc 

L,M,N,0 

M 


Me 

Mz 

Mq 


n 

P 

q. 


s,t 


U 

x,y 

3 


functional determinant, or Jacobian 

transonic similarity parameter 

critical value of similarity parameter 

potential limit of similarity parameter 

constants, identified by superscripts and subscripts 

local Mach number 

critical Mach number 

potential limit Mach number 

free-stream Mach number 

number of approximations 

particular integral 

function of s 

elliptic coordinates 

free-stream velocity 

cartesian coordinates 

'v/l - 
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r„ Mo2(i + Vl 

y ratio of specific heats 

A Laplacian 

6 constant 

1 ] transformed y-coordinate 

T thickness ratio 

<P peirturhation potential 

consteints, identified hy suhscrlpts 

Subscripts ; 

indicate partial derivatives 
h, J,k,m,p,r,v summation indices 
Superscripts ; 

X,n,v number of approximation 
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ZASIZ 17 


Klltfpg APywIT WtTTl'W yCB TT.T.TPfig 


4 

M A 


hj 


b \ 

2 

5 

.4 

5 ^ 

6 

7 

6 

9 

10 

1 

0.14320 

0.09286 

- 0.16501 

0.04965 

0.01266 

0.00436 

6.00179 



0.00007 

3 

- 0.21460 

- 0.43555 

, 0.31663 

0.02354 

- 0.02353 

- 0.02820 

- 0.00936 

— 

- 0 . 00007 ' 

5 

0.07160 

0.57662 

0.19904 

- 0.52946 

0.06649 

0.02305 

0.03632 

— 

0.00024 

7 

— 

- 0.26551 

- 0.75659 

0.80076 

- 0.07486 

0.05956 

- 0.07065 

— 

- 0.00132 

9 

— 

0.02977 

0.49494 

- 0.38547 

0.02825 

- 0.10221 

0.06084 

— 

0.00262 

11 

— 

— 

- 0.09795 

0.04100 

— 

0.04608 

- 0.01895 

— 

- 0.00222 

13 

— 

— 

0.00693 

— 

— 

- 0.00464 

— 

— 

0.00069 





4 

1 

0.16250 

3 

- 0.48749 

5 

0.48749 

7 

- 0.16250 



1 

2 

5 

4 

5 

6 

7 

B 

9 

10 

11 

12 

13 

1 

0.22496 

- 0.09547 

- 0.02332 

o.oioia 

- 0.00414 

- 0.00072 

- O.OOOIB 

- 0,00006 

— 

- 0.00000 

— 

— 

— 

3 

— 

- 0.08592 

- 0,04018 

0.06871 

- 0.00294 

0.00241 

0.00141 

0.00034 

— 

0.00000 

— 

— 

— 

S 

~ 

0.00682 

0,07205 

- 0.05206 

0.01638 

- O . OL 209 

- 0.00092 

- 0.00188 

— 

- 0,00001 

“ 

“ 


7 

— 

— 

- 0.01328 

- 0.04347 

0.00673 

- 0.01152 

- 0.01829 

0.00942 

— 

0.00005 

— 



9 

— 

— 

0.00085 

0.01523 

- 0.01377 

0.00125 

- 0.00639 

0.00716 

, - 0.00132 

- 0.00028 

— 

— 


U 

— 

— 

~ 

- 0.00187 

0.00085 

— 

0.00134 

- 0.00066 

~ 

- 0.00021 

0.00012 

— 

— 

15 

— 

— 

— 

0.00009 

— 

— 

- 0.00008 

— 

— 

0.00002 

— 

— 

- 0.00000 


(4)Mjg 



3 

4 

5 

6 

7 

1 

. 

- 0.02167 

— 

— 

— 

3 

- 0.14445 

0.13928 







5 

— 

0.10833 

- 0.10590 

— 

— 

7 

— 

- 0,00985 

— 

— 

0.00851 
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mATJTTP XT 


(a) 




2 

3 

4 

5 

A 

7 

3 

9 

10 

11 

T O 

13 

1 j 
J.V 

1 

-0.30193 

-0.19052 

-0.04166 

0.16206 

0.01136 

0.00283 

0.00531 

0.00204 

0.00196 

0.00004 

0.00028 



3 

0.45289 

0.88909 

0.58486 

-0.74174 

-0.U23B 

-0.02831 

0.02498 

-0.05838 

-0.00891 

0.00041 

-0.00066 

___ 



5 

-0.1 FWW 

-1 .IflSSI 

-1 .Sfinsfi 

1.18749 

0. 75422 

0.01647 

-0.27262 

0.16464 

0.02092 

0.00200 

0.00086 





7 


0.54470 

2.49445 

0.21268 

-2.59319 

0.34671 

0.49167 

-0.15302 

-0.05639 

-0.01005 

-0.00175 




9 



-0.06107 

-1.32017 

-2.07283 

3.52376 

-0,63624 

-0.11714 

-0.21780 

0.10470 

0.00985 

0.00719 

- 

0.00003 

ai 





0.28125 

1.67598 

-1.93465 

0.33485 

-0.34649 

0.43170 

-0.09089 

0.00799 

-0.01371 


-0.00015 

13 

___ 

— 

-0.01847 

-0.48297 

0.57619 

-0.03630 

0.25889 

-0.20890 

0.02862 

-0.01837 

0.01153 

— 

0.0002,6 

15 

— 

— 


0.06266 

"0*02552 


-0.04744 

0.01972 

— 

0.00892 

-0.00364 

— 

-0,00019 

'17 

— 

— 

— 

-0.00327 

— 

— 

0,00290 

— 

— 

-0.00079 

— 

— 

0,00006 


' (Wljj 




4 

5 

6 

7 

8 

0 

m 

0.45334 

0.14341 

0.19737 


0.00516 


1.30002 

0.06802 

-1.06659 



5 

-1.50002 

-1.52956 

2.24554 

— 

0.17670 

■ 

r\ 1 vwA 


O rti 

— 

f\ r 

“V*\T*XVO 

9 




bbi 

0.28668 

U 

— 


— 

— 

-0.08735 
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TABM T - ITETH JPPEOnMATIO¥ TOE mjPBB - Oattolnded 

(o) 


iiX 

1 

2 

3 

4 

5 

E 

7 

e 

e 

10 

u 

12 

IS 

14 

16 

17 

1 

0.2€6fi0 

-O.lOOSl 

-0.02381 

-0.01818 

0.00824 

0.00228 

0,00006 

0.00010 

0.00005 

0.00002 

O.OOOOO 

0,00000 

— 

— 

— 

— 

3 

... 

-0.09058 

0.31449 

0.29580 

-0.04370 

-0.0E153 

-0.00071 

O.OOOE4 

-0.00063 

-o.oooao 

0.00000 

-0,00000 

-u- 

, 

... 

— 

6 

... 

0.0071B 

0.07389 

0.08941 

-0.41898 

0.29128 

0,00069 

-0.00(513 

0.00294 

0.00028 

0,00002 

0.00000 

— 

— 

... 

— 

7 

... 

... 

-0.01392 

-0.07241 

0,03120 

0.05421 

-0,OC332 

0.00882 

-0.00478 

-0.00111 

-0.00014 

-0.00002 

— 

— 

— 

— - 

9 


— 

o.oooas 

0.02051 

0.03343 

-0.07440 

0,0193G 

0.02276 

-0.020S2 

0.00551 

0,00025 

O.OOOLL 

— 

0.00000 

— 

— 

u 


— 

... 

-0.00249 

-0,01733 

0.03878 

-0.00485 

0.00566 

-0.01079 

0.00453 

0.00076 

-0,00060 

— 

-0.00000 

— 

... 

13 

— 

— 

— 

O.OOOIE 

0.00535 

-0.00223 

0,00050 

-0.00248 

0.00237 

-0.00042 

0.00058 

-0.00046 

0.00010 

0.00001 

— 

— 

is ’ 

— 

— 

... 

... 

-0.00031 

0.00015 

— 

0.00030 

-0.00014 

— 

-0,00009 

0.00004 


0.00001 

- 0.00000 

— 

17 

— 

— 

... 

— 

0.00001 

— 

— 

- 0,00001 

— 

— 

0,00000 

— 

— 

- 0.00000 

— 

- 0.00000 


(a) «hj 



3 

4 

5 

6 

7 

8 

9 

10 

n 

1 



-0.00689 

0.00698 

0.00584 

— 

o.ooooe 

— 

— 


3 

-0.14818 

0.18578 

0.00126 

-0.03613 

... 

-0,00073 

... 

— 

— 

8 


0.14446 

-0.1031S 

0.20414 

--- 

0,00453 

--- 

— 

— 

7 

— 

-0.01313 

-0.09638 

0,15737 

-0.02476 

-0.02274 

— 

— 

— 

B 

— 


0,01988 

— 

— 

-0.01688 

0.01210 


— 

11 

— 

— 

-0.00123 

— 

— 

0.00154 

— 

... 

-0,00038 



6 

(») IT hj 
H ®55 - 0,076ia 
N®j^ - 0| h,J / S 
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TABLE VI - SIXTH APPEOHMAnON PC® ELLIPSE 

(a) A 


h \ 

2 

3 

4 

5 

6 

1 

7 

8 

9 

1 

- 0.35512 

- 0.20367 

- 0.54907 

0.53760 

0.52224 

0.39212 

0.01983 

- 0.00962 

3 

0.53268 

0.95047 

2.16610 

- 1.13329 

- 3.53878 

2.24494 

0.27821 

- 0.05669 

5 

- 0.17756 

- 1.26381 

- 3.67930 

- 0.80210 

8.32920 

- 3.82950 

- 0.96660 

0.19693 

7 



0.58229 

3.23792 

4.50042 

- 8.29361 

0.43150 

1.86413 

0.21831 

9 



- 0.06528 

- 1.43991 

- 5.29265 

1.98519 

6.03432 

- 2.58489 

- 0.95363 

11 




0.28437 

2.85636 

2.82542 

- 7.73804 

2.34101 

0.50963 

13 

— 

— 

- 0.02010 

- 0.75951 

- 2.64252 

4.22655 

- 1.11805 

0.50443 

15 

o 


— 

0.09828 

0.97957 

- 1.11033 

0.22088 

- 0.54201 

17 




- 0.00511 

- 0.18374 

0.13964 

- 0.01486 

• 0 . 14955 ' 

19 

— - 

— 

— 

— - 

0.01776 

- 0.00697 


- 0.01773 

21 

— 

““ 

— 

— 

- 0.00072 

... 


0.00082 


6 

(a) A jjj - Caijcluded 


n \ 

10 

11 

12 

13 

14 

15 

16 

18 

1 

0.00467 

- 0.00511 

0.00088 

- 0.00052 

- 0.00014 

- 0.00002 

- 0.00000 

- 0.00000 

3 

0.00801 

0.01975 

0.00047 

0.00045 

0.00043 

0.00006 

0.00002 

0.00000 

5 

- 0.02668 

- 0.05440 

- 0.01173 

0.00009 

- 0.00080 

- 0.00001 

- 0.00007 

- 0.00000 

7 

- 0.33773 

0.17216 

0.01806 

0.00866 

0.00131 

- 0.00007 

0.00008 

0.00000 

9 

0.69461 

- 0.17444 

0.00844 

- 0.03344 

- 0.00290 

- 0.00018 

- 0.00005 

- 0.00000 

11 

0.11262 

- 0.17842 

- 0.02885 

0.03281 

0.00916 

0.00084 

0.00007 

0.00000 

13 

- 1.00872 

0.39807 

- 0.01228 

0.02266 

- 0.01671 

- 0.00085 

- 0.00029 

0.00001 

15 

0.66684 

- 0.19573 

0.05158 

- 0.05514 

0.01389 

- 0.00039 

0.00059 

- 0.00003 

17 

- 0.12088 

0.01812 

- 0.03191 

0.02680 

- 0 .. 00423 

0.00115 

- 0.00053 

0.00004 

19 

0.00725 


0.00568 

- 0.00239 



- 0.00057 

0.00017 

- 0.00003 

21 

— 


- 0.00033 



— 

0.00005 


0.00001 
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TABIS 71 - SIXTH APPEOHMATIC® FOE ELLIPSE - Continued 

0^) B hj 



4 

5 

6 

7 

'8 

9 

10 

12 

1 

-0.44175 

0.27651 

0.17244 

-0.21716 

-0.07427 

0.00325 

-0.01077 

-0.00071 

3 

1.32524 

0.07448 

-0.95724 

1.36008 

0.22528 

0.09495 

0.02864 

0.00142 

5 

-1.32524 

-2.37617 

1.17188 

-1.02543 

-0.60779 

-0.40598 

-0.11295 

-0.00071 

7 

0.44175 

3.51929 

1.55575 

-3.43737 

1.2i433 

0.34626 

0.45626 

o 0.00399 

9 



-1.66999 

-3.93061 

5.82295 

-1.13510 

0.47585 

-0.84175 

-0.02074 

11 


0.17588 

2.42446 

t2. 77921 

0.37754 

-0.89299 

0.69364 

0.02897 

13 



— 

-0,46759 

0.27615 



0.41725 

-0.20907 

-0.03169 

15 

— 

— 

, 0.03092 

- — 

— 

-0.03859 

— 

0.00947 





1 

-0.38239 

3 

1.24277 

5 

-1.43397 

7 

0.66918 

9 

-0.09560 


f 


2275 





2275 


5 


MCA aasf 2511 


33 




TflELE 71 - SnXE APHROnMAnCffl TOR ELUPSE - Contlmed 


(d) thj 


\ J 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1 

0.34464 

- 0.11837 

- 0.02546 

- 0.05231 

0.02720 

0.01467 

- 0.00949 

- 0.00061 

- 0.00011 

0.00003 

3 



- 0.10654 

- 0.40107 

0.52581 

- 0.06792 

- 0.13205 

0.06802 

0.00625 

- 0.00046 

0.00016 

5 



0.00846 

0.07899 

0.27792 

- 0.66704 

0.58684 

- 0.18449 

- 0.03118 

0.00119 

- 0.00076 

7 

— 

— 

- 0.01456 

- 0.09487 

- 0.12642 

0.65280 

- 0.67976 

0.21063 

0.01291 

- 0.00314 

9 





0.00091 

0.02215 

0.08919 

- 0.06701 

- 0.10896 

0.08921 

0.00838 

- 0.01008 

11 







- 0.00271 

- 0.02956 

- 0.02921 

0.09997 

- 0.03921 

- 0.02279 

0.02609 

13 







0.00013 

0.00527 

0.01955 

- 0.03495 

0.01065 

- 0.00476 

0.01374 

15 



... 

... 



- 0.00049 

- 0.00517 

0.00631 

- 0.00137 

0.00373 

- 0.00533 

17 

... 





... 

0.00002 

0.00073 

- 0.00058 

0.00007 

- 0.00064 

0.00063 

19 

... 

... 


... 

... 

- 0.00005 

0.00002 


0.00006 

- 0.00005 

21 

— 


— 

... 

... 

0.00000 

... 

_ „ _ 

- 0.00000 

... 


(d) L jjj - Ccmoluded 



U 

12 

•13 

14 

15 

16 

17 

18 

19 

21 

1 

- 0.00004 

0.00001 

- 0.00000 

- 0.00000 

- 0.00000 

- 0.00000 

___ 

- 0.00000 



3 

0.00018 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

... 

0.00000 



5 

- 0.00057 

- 0.00010 

0.00000 

- 0.00000 

- 0.00000 

- 0.00000 


- 0.00000 



7 

0.00239 

0.00020 

0.00007 

0.00001 

- 0.00000 

0.00000 

___ 

0.00000 



9 

- 0.00436 

0.00001 

- 0.00038 

- 0.00003 

- 0.00000 

- 0.00000 



- 0.00000 

... 


11 

- 0.01128 

0.00051 

0.00068 

0.00012 

0.00001 

0.00000 


0.00000 



13 

- 0.00829 

- 0.00073 

0.00227 

- 0.00062 

- 0.00002 

- 0.00000 

... 

0.00000 



15 

0.00188 

- 0.00060 

0.00098 

- 0.00048 

0.00001 

0.00002 

... 

- 0.00000 



17 

-O.OOOU 

0.00022 

- 0.00022 

0.00005 

- 0.00002 

0.00002 

- 0.00000 

0.00000 



19 



- 0.00003 

0.00001 



0.00000 

- 0.00000 



o.ooboo 

0.00001 


21 

— 

0.00000 

— 

— 

- 0.00000 

— 

— 

- 0.00000 



0.00000 



01 

!>■ 


aiABLa 71 - anim APPEonnAiniQir toe ttt.t.tpsh: - coaoiadod 


6 

(e) M 





3 

4 

6 

6 

7 

3 

9 

10 

11 

12 

13 

15 

1 


- 0,02945 

0.01152 

- 0.00256 

- 0.00452 

- 0.00118 

0.00004 

- O.OOOll 



- 0.00000 

— 


3 

- 0.15841 

0 , ia 032 

0.00465 

0.00546 

0.05400 

0.00410 

0.00132 

0.00031 

— 

0.00001 

— 

— 

5 


0.14725 

0.10397 

- 0.17789 

- 0.04273 

- 0.01558 

- 0.00725 

- 0.00151 

— 

- 0.00001 

— 

— 

7 

«« p . 

- 0.01339 

- 0,14664 

- 0.11967 

- 0.02132 

0.04857 

0.01062 

0.00887 

— 

0.00004 


— 

9 




0.02982 

0.06621 

- 0.18197 

0,06677 

0,05151 

- 0.04430 

— 

- 0.00053 

— 

— 

11 



- 0.00183 

- 0.02852 

0.05860 

- 0.00662 

0,02232 

- 0.03305 

0.00811 

0.00169 

— 

— 

13 


_ 


0.00362 

- 0.00230 

— 

- 0.00474 

0.00303 

— 

0.00127 

- 0.00087 


15 

— 

— 

— 

- 0.00016 

— - 

— 

0.00027 


— 

- 0.00012 

— 

0,00001 



\-5 

h\ 

5 

6 

7 

9 

1 

— 

- 0.01093 

— 

— 

3 

— 

0,04603 

— 

— 

5 

0.11881 

- 0.15036 

— 

— 

7 

— 

- 0.05148 

0.12278 

— 

9 

— 

0.00212 

— 

- 0.01322 
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TABLE 711 - FIRST APPBOZIMATIOU FOE KAPLAE SKCTICm 

1 

(s) 


(b) 


\ J 
b\ 

1 

3 ■ 

1 

0.75000 

— — 

3 

— 

-0.25000 


TABLE TUI - SECOHD APPEOXIMATION FOR KAPLAN SECTION 


<=> * w {if 



2 

4 

6 

1 

-2 

-3 

-1 

3 

3 

' 9 

2 

5 

-1 

-9 

-4 

7 

0 

3 

5 

9 

0 

1 

0 

-2 


(b) L 


bj 


e-s 

1 

2 

3 

4 

5 

6 

7 

9 

1 

0.32412 

-0.09375 

— 

-0.02813 

— 

-0.00402 

. 

— 

3 

— 

-0.08438 

-0.20566 

0.18081 

— 

0.01042 

— 

— 

5 

— 

0.00670 

— 

0.14063 

-0.05382 

-0.05114 

■ 

— 

7 

— 

— 

— 

-0.01278 

— 

-0.05409 

0.05367 

— 

9 

— 

— 

— 

— 

— 

0.00625 

— 

-0.00417 





lABii EC - TEn® APEFormAnoa job kaplw BBjnoa 


OI 

<J) 


(«) A ijj 


ll\ 

£ 

s 

4 

B 

6 

7 

6 

a 

10 

11 

12 

I 

-0,2i306 

-0.07BIB 

-0.34704 

-0.07433 

-O.lfiOfil 

-0.26697 

-0.20181 

0.06130 

-0.14087 

___ 



I 

0.56462 

0.56565 

1.22544 

0.14304 

-0.0061E 

0.43360 

-0.49405 

0.24415 

-0.55548 

--- 

-0.02109 

s 

-0.12154 

-0.46606 

-1.54496 

-1.41706 

0.17637 

-0.55128 

1,02705 

“0,35314 

0.27S4S 

... 

0.07754 

7 

— 

o.esioo 

0.46358 

2.10642 

0.16590 

1.43764 

-1.80276 

0,20793 

-0.59863 

... 

-0.09141 

9 

— 

-0.02C10 

— 

-1.0177B 

-0.14624 

-2.49763 

1.70192 

-0.37881 

0.96463 


0.05526 

U 

— 

— 

--- 

0.1090a 

--- 

L. 32190 

-0.S5596 

0.87676 

-1.06352 

— 

-0.096U 

13 

— 

_ 

— 

— 


-0.14718 

— 

-0.EE56B 

-0,13069 

... 

0.1I2S0 

IS 

— 

— 

--- 

— 

--- 

— 

— 

0.06326 

— 

--- 

-0.04219 


(b) L^j 


XJ 

ll\ 

1 

s 

‘3 

4 

6 

6 

7 

8 

8 

10 

11 

12 

13 

14 

1 

0.39487 

-0.06103 

-0.00979 

-0.0S514 

0,00310 

-0.00617 

-0.0060E 

0.00320 

0.00077 

0.00142 



mmm 

MM 

3 

— 

-0.(77282 

-0.21244 

0.17476 

0,00694 

-0.00023 

0.00215 

-0.00098 

0.00339 

0.00619 

--- 

-0,00016 

mmm 



5 

— 

0.00679 

0.03039 

0.149U 

-0.15715 

0.01622 

-0,00500 

0.02633 

-0,00938 

0.00967 



0.00065 

mmm 



7 

— 

— 

-0.00360 

-0.01420 

-0.08735 

-0.01184 

0.22204 

-0.12S52 

0.00650 

-0.01174 

..... 

-0,00096 

... 



fi 

— 


O.OOOK 


0,01617 

0.00525 

0.07006 

-0,10011 

-0.03949 

0,06077 


0.00100 

... 




— 

— 

— 

— 

-0.00114 

— 

-0.00636 

0,00972 

-0.02197 

0.05063 

-o.oieaa 

“0.00428 

... 



13 

— 

— 

— 

— 

— 

— 

0.00123 

— 

0.00595 

O.OOIBS 

... 

-6.00460 

-o.ooEoe 

— 

15 

— 

— 

— 

— 

— 

— 

— 

— 

-O.OOOU 



... 

0.00052 

... 

-0.00015 




s 

3 

s 

7 

9 

3 

-0,06094 






5 



0.14170 




7 



■ ,,,, 

-0.10415 

— 

_fl 

T, 

- 

... 

Q-oaios 


QiSB 
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TABia I - rOaHTH AEPROTIMAIIOT TOH KAPLAIT aKGEEOn 


4 


\ J 

2 

3 

4 

5 

6 

7 

8 

9 

1 

-0.27369 

-0,10158 

-0.70184 

0.1S38i 

-0.25999 

-0.08941 

0.27562 

0.00827 

3 

0.41053 

0.47402 

2.21465 

0.20698 

-0.14470 

0.34639 

-0.94946 

0.74600 

5 

-0.13684 

-0.63029 

-2.54101 

-2.37393 

-0.09017 

-0.37686 

3.50541 

-2.30284 

7 

— 

0.290^ 

1,28382 

3.50709 

2.64389 

-6.52365 

1.77718 

-2.32167 

9 

— 

-0.03256 

-0.31320 

-1.66637 

-3.74815 

-3.76560 

2.16703 

1.62381 

11 

— 

— 

0.06195 

0.18240 

1,95475 

2.08185 

4.82315 

-3.15023 

13 

— 

— 

-0.00437 

— 

-0.38139 

-0.24116 

-4.05218 

1.40343 

15 

— 

--- 




• 0.02576 

— 

0.80760 

-0.10562 

17 

— 

— 

— 

— 

— 

— 

-0.05443 

— 

19 

— 

— 

— 

— 

— 

— 

— 

— 

21 

— 

— 

— 


— 

— 

— 

— 


4 

(a) A jjj - CcExolnded 


h\ 

10 

11 

12 

13 

14 

15 

16 

18 

1 

0.30919 

-0.13252 

-0.03357 

-0,00319 

-0.03780 

0.003U 

-0.00705 

-0.00035 

3 

-1.56053 

0.63407 

-0.17619 

0.50166 

0.04789 

0.00980 

0.01409 

0.00035 

5 

1.92675 

-0.67733 

0.58936 

-0.64026 

-0.13063 

-0.06476 

-0.00705 

... 

7 

1,37105 

1.37105 

-0.18288 

0.12860 

0.42178 

0.09728 

-0.02032 

... 

9 

5.14985 

-2.08886 

-0.71016 

-0.22855 

-0.46413 

0.00347 

0.071U 

-0.00280 

11 

-2.88950 

-2.33567 

0.76885 

1,04169 

0.41837 

-0.03373 

-0.07825 

0.00796 

13 

-3,37322 

6.15880 

-0.99893 

0.85868 

-0,95828 

-0.07389 

0.02888 

-0.00818 

15 

3.42262 

-3.15298 

1.55071 

-2.77829 

1.12471 

-0.09933 

-0.01642 

0.00786 

17 

-0.71285 

0.22345 

-1.01520 

1,34638 

-0.42191 

0.32418 

0.02302 

-0.01576 

19 

0.04935 

— 

0,22411 

-0.02672 

— 

-0.18607 

-0.00802 

0.01698 

21 

— 

— 

-0.01630 

— 

— 

0.01994 


-0.00606 




ll\ 

4 

6 

8 

10 

12 

1 

0.01714 

0.02471, 

0.02657 

0.00976 



3 ' 

-0.05142 

-0.11585 

-0.07267 

-0,01081 

-0.00395 

5 

0.05142 

0.24498 

0.17382 

0.07426 

0.01578 

7 

-0.01714 

-0.22984 

-0.35115 

-0.09793 

-0.01973 

9 

— 

0.07599 

0,33867 

0,18995 

0.01776 

u 

— 

— 

-0.11524 

-0.19764 

-0.03354 

IS 

— 

— 



0,07041 

0.03749 

15 

— 

— 

— 

— 

-0.01381 
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TABI£ X TC Smm AFFBtUQKASIQZf 7QH EAPLAU SECTION - Cosiolodod 

4 

- (o) hj 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1 

0.37947 

-0.09123 

-0.01270 

-0.04618 

0.00685 

-0.00719 

-0.00186 

0.00450 

0.00010 

0.00314 

3 

— 

-0.082H 

-0.36592 

0.30799 

0.01294 

-0.00727 

0.00866 

-0.01764 

0.01036 

-0.01727 

5 

— 

0.00652 

0.05959 

0.28741 

-0.30615 

0.01610 

-0.01570 

0.09171 

-0.04112 

0.02595 

7 

— 

' 

-0.00726 

-0.03905 

-0.14613 

-0.21970 

-0.78172 

-0.45988 

0.05554 

-0,04627 

9 


— 

0.00045 

0.00482 

0.03011 

0.08574 

0.11761 

-0.10872 

-0.30580 

0.27104 

n 





— 

-0.00059 

-0.00190 

-0.02300 

-0.02891 

-0.08519 

0.07876 

0.12864 

13 

— 

— 



0.00003 

— 

0.00287 

0.00201 

0.05859 

-0.0159S 

0.04919 

IS 

— 

— 

— 

— 


-0.00014 

— 

-0.00502 

0.00073 

-0.02738 

17 

— 

— 

— 

— 



— 

— 

0.00024 

— 

0.00377 

19 





— 













-0.00019 

21 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 


4 

(c) L - Conolaflftd 


s 

11 

12 

13 

14 

15 

16 

17 

18 

19 

21 

1 

-o.oono 

-0.0Q023 

-0.00002 

-0.00019 

0.00001 

-0.00005 

— 

-0.00000 

— 

— 

3 

0.00566 

-0,00150 

0.00169 

0.00026 

0.00005 

0.00006 

— 

0.00000 

— 

“ 

5 

-0.00706 

0.00498 

-0.0044S 

-0.00076 

-0.00082 

-O.OOOU 

— 

— 

— 

— 

7 

0.01904 

-0.00198 

0.00107 

0.00287 

0.00035 

-0.00010 

— 

— 

— 

— 

9 

-0.0S222 

-0.01127 

-0.002SO 

-0.00404 

0.00002 

0.00041 

— 

-0.00001 

— 

— 

11 

-0.14319 

0.05191 

0.02170 

0.00558 

-0.00052 

-0.00058 

— 

0.00004 

— 

— 

IS 

-0.12831 

0.04140 

0.05407 

-0.05549 

-0.00132 

0.00033 

— 

-0.00005 

— 

— 

IS 

0,05032 

-0,01919 

0.04961 

-0.05878 

0.00756 

-O.OOOS3 



0,00008 



— 

17 

-0.00133 

0.00700 

-0.01122 

0.00454 

-0-00507 

-0.00070 

0.00404 

-0.00045 





19 



-0,00105 

O.OOOU 



0.00157 

0.00008 



-0.00046 

-0.00006 



21 

— 

0.00008 

— 

— 

-0.00009 

— 

— 

0,00005 

— 

-0.00 002 


(4) « hj 



3 

4 

5 

6 

7 

6 

9 

10 

u 

12 

13 

IS 

1 

. — 

O.OOU4 

— 

0,00071 

— 

0,00042 



0.00010 

— 

— 

— 

— 

3 

-0.07900 

-0.00735 

— 

-0.00429 

— 

-0.00132 



-0.00054 

— 

-0.00003 





5 
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Tlgnre 2. - ParsboUo ourre of ■transonic eq.ua'tlon. 
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Figure 4. - Functional determinant of transformation 
I = cosh s cos t, Tj = sinh s sin t. 



1 - 

l_. 


Figaro S. - Kaplan soctlon ani region of -yalldlty of olmpUflod flow equation. Thlokneaa ratio, O.lj froo-streain 

Maoh nisaber,- 0*9. 
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Figure 8. - Estimation of potential limit similarity parameter for 

elliptic cylinder. 
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Figure 9. - Maximum velocity increment on surface of thin ellipse. 



Figure 10. - M a ximum local Mach number m surface of thin ellipse. Six 

approximations . 
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(a) Similarity parameter, 0.3. 



(b) Similarity parameter, 0.4. 
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(o) similarity peurameter, 0.5. 

Figure 11. - Telocity Increment at surface of ellipse. 
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Figure 16. - Maximum local Mach number on aurfaco of Kaplan section. 
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(a) Similarity parameter, 0.2. 

Figure 17. - Velocity increment at surface of Kaplan section. 
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I = cos t 

(c) similarity parameter, 0.4. 

Figure 17. - Concluded. Telocity increment at surface of Zaplan section. 
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